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Abstract 

Based on the revealed hidden supergroup structure, we develop a new approach 
to the Hubbard model. We reveal a relation of even Hubbard operators to the spinor 
representation of the generators of the rotation group of four-dimensional spaces. We 
propose a procedure for constructing a matrix representation of translation generators, 
yielding a curved space on which dynamic superfields are defined. We construct a new 
deformed nonlinear superalgebra for the regime of spinless Hubbard model fermions 
in the case of large on-site repulsion and evaluate the effective functional for spinless 
fermions. 
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1 Introduction 



The Hubbard model [1] remains the main test ground for investigating the effects of strong 
correlations between electrons. Traditionally, the effects of a strong Coulomb interaction play 
an important role in understanding the mechanisms of high-temperature superconductivity, 
the physics of the Mott-Hubbard-type metal-insulator transition, and the related magnetic 
states. Experimental evidence of a spin liquid in organic metal (BEDT-TTF)2Cu2(CN)3 [2] 
additionally complicated the situation: it has become necessary to take the spin liquid into 
account in this model. We believe that this has revitalized the problem of developing new 
approaches to the Hubbard model based on field theory methods and especially on using the 
functional integral techniques. 

Using the path integral proposed for the Hubbard model in [3], [4], we here continue 
developing the superfield formulation. We discover a nonlinear hidden supergroup occurring 
in strongly correlated systems and show how a matrix representation for the coordinate shift 
generators can be introduced together with the collective dynamic fields. We compute the 
effective functional based on the new nonlinear superalgebra, also obtained here. 

We recall the crucial points of the approach in [3]. As the starting formulation, we take 
the atomic representation of the Hubbard model. A property of this description is that it 
allows introducing a local supergroup whose generators arc defined by the full set of the on- 
site Hubbard operators. The supergroup generators act globally and are independent of the 
lattice coordinates. Next, a supercoherent state is constructed that specifies a superorbit of 
our supergroup and contains a set of dynamic (i.e., time and coordinate dependent) fermionic 
and bosonic fields. Passing from the operator to the superfield formulation is achieved using 
the effective functional of the form 

f _ {G\d/dT - ffHub|g) ... 

{G\G) ' 

where \G) is the supercoherent state that specifies the supergroup orbit to be defined below 
and i^Hub is an operator expression for the Hubbard model. 

We emphasize that one of the fundamental questions of strongly correlated systems is 
the problem of the existence of a hidden supergroup structure and the procedure for intro- 
ducing it into the Hubbard model. Identifying this superstructure would allow clarifying the 
question of new symmetries and supersymmetries and developing a procedure for calculat- 
ing effective superfield functional for the strong Coulomb interaction systems described by 
models that are nonlinear in the dynamic bosonic and fermionic fields. 

A central element of the proposed approach is the procedure for identifying even Hubbard 
operators with local generators of the Lorentz group and other rotation groups of four- 
dimensional spaces. In strongly correlated systems, we understand the Lorentz group to be 
a four-dimensional rotation groTip containing hyperbolic rotations with a "speed of light" 
constant depending on the characteristic energy, for example, the parameter of the electron 
band width or the exchange integral. The full set of Hubbard operators is identified with a 
Lorentz group superextension, which is a subgroup in the superconformal group. 

Another important point relates to the introduction of a matrix representation for the 
group of translations. This subgroup, via a quotient, allows introducing the coordinate space 
on which the dynamic bosonic and fermionic fields are defined. We introduce the subgroup 
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of translations using the procedure not of the group affinization, as in the majority of papers 
(see, e.g., [5]), but of the deformation and contraction of algebras. Our approach is more 
complicated but seems to be the only one suitable for the Hubbard model. The tricks that we 
use allow evaluating the effective functional for some deformation of the algebra of spinless 
fermions, which is nonlinear in the Fermi and Bose generators. 

We show how the hidden supergroup structure emerges in the Hubbard model. Any 
quantum system is characterized by the wave function ^, which is a function of the coordi- 
nates X, y, and z and time t, i.e., = vl/(a;, z,t). The most widespread groups act on this 
wave function via operators defined in differential form. For example, one of the widespread 
continuous groups acting on the wave function is the Poincare group, which contains the 
space-time translation operators 

. d . d . d Id 

OX oy oz c at 

the operators of spatial rotations in the {x\ x^) planes 

d d\ , f d d\ , f d d 
ox oy J \ oz ox J \ oy oz 

and the operators of hyperbohc rotations (boosts) in the(ct, x'') planes 



f d d\ , f d d\ , f d d 

The action of a finite group on the wave function is given by the exponential map and has 
the form 

gBiPi+^f.-L.+i/fcEfc^^^^ ^ e-^'=^'=*(Ax + B), (2) 

where Ei and Hk are small parameters, summation over the indices A; = 1, 2, 3 is under- 
stood, the matrix A defines a vector representation of the rotation group, and the matrix B 
represents shifts. We note that the parameters determine both the matrix A acting on 
the coordinates and the elements of the rotation group for the function ^ in the spinor basis. 
The wave function is defined in the spinor representation and is acted upon by elements 
of the group parameterized by the H^. 



2 Hubbard model 

The Hubbard model [1] , [6] is one of the fundamental models describing systems of strongly 
interacting electrons in solids. It underlies the description of the band magnetism and 
superconductivity in strongly correlated metals and the metal-insulator phase transition in 
sohds. 

The Hubbard model involves electrons that are nondegenerate with respect to the orbital 

state and propagate by hops over sites of the crystal lattice with a Coulomb interaction at 
one site. In the secondary quantization representation, the Hubbard model Hamiltonian has 
the form 

H = -W al^ttja + rii^nn - fi^riia, (3) 

i,j,(7 i i,a 
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where and Ojo- are Fermi operators of creation and annihilation of an electron with the 
spin a at the ith site, — a^^Qiu is the operator of the number of electrons with a given 
spin a at the ith site, and the spin a takes two values t (+) and I {—). 

The Hubbard model involves only three parameters: the matrix element W of the electron 
transition to the neighboring site, which is related to the electron band width, the parameter 
U of Coulomb repulsion of two electrons at a site i, and the chemical potential // or the 
electron concentration n (the average number of electrons per lattice site). 

In what follows, we set W = 1 and consider the half-filled band in the case of strong 
repulsion. In that regime, the lowest-lying states are described by magnetic fields together 
with fermionic degrees of freedom. The repulsion parameter practically drops out of the 
problem. 

We are interested in the two-dimensional Hubbard model, and we start investigating it 
from the atomic limit. 



3 Hubbard operators 



Under the condition of a strong Coulomb repulsion, the Coulomb term in the Hamiltonian 
is taken as the zeroth approximation. Then the zeroth approximation reduces to a one-site 
problem and can be solved in the basis of localized atomic functions: 



\ip): |iO), 



(4) 



where \iO) is the state in which the ith site does not contain an electron, \i+) = is the 
state with one spin-up electron at the ith site, |i— ) = is the state with one spin-down 
electron at the ith site, and |i2) = |iti) is the state with two electrons at the ith site, one 
with spin up and the other with spin down. Any state at the site i can be represented as a 
superposition of these localized functions: 



ipi = ai\iO) + /3i\i+) + 7i|i-) + Si\i2), 



(5) 



where ctj, 7^, and are coefficients of the atomic basis functions with /5j and 7^ being 
Grassmann numbers. In what follows, we use matrix elements of the superspinor represen- 
tation of the local supergroup to represent the coefficients of this decomposition. 

Transitions between the atomic basis states are described by 4x4 matrices corresponding 
to the Hubbard operators: 



X: 



pq 



\ip){iq\ 





xr 


xr 


xr\ 


xr 


xr 


xr 


xf 


xr' 


x-^ 


X-- 


xf 


\x!' 




xr 


xf) 



(6) 



Each Hubbard operator corresponds to a matrix with the unity at the intersection of the pth 
row and the gth column and with all other elements equal to zero. The Hubbard operators 
satisfy the relations Xf^X['^ = bqrX\^ ^ which define their algebra, and the completeness 
condition = 4 for the diagonal elements. 
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The Fermi (f-type) operators are X2+, and X^" because 

they change the number of electrons at a site by an odd number, 1 or 3. They satisfy the 
anticommutation relations 

[Xf , XY]+ = + = Sij {s,rxr + SgrXl'^). (7) 

The Bose (b-type) operators are , X~~^ , Xf^, and Xf^. They change the number of 
electrons at a site by an even number, or 2. They satisfy the commutation relations 

[Xf , Xf ]_ = Xf^Xf - Xf Xf = 5i,(V^r - (8) 

This formula also holds if it involves at least one b-type operator or one of the diagonal 
operators X™, X++, X-— , or Xf . 

In terms of the X operator, the Hubbard model Hamiltonian has the form 

H = 5^{-/.X++ - ^Xr- + {U- 2/x)Xf } - 

i 

- W 5^{(Xr + X^)(X;° + X-2) + (Xro + Xf^){X^- + X;^)}, (9) 

ij 

where the first sum also inchidcs the chemical potential fi. In this representation, the 
Coulomb term is linear, and the kinetic energy term is given by a bilinear form in the 
X operators. In what follows, we use this operator formulation because it takes us precisely 
to the supergroup construction, which, in our opinion, is the central element of all strongly 
correlated models. 



4 Minimal superalgebra in the Hubbard model 

We use the construction of a "symmetry tower" introduced for the Hubbard model in [4]. 
We consider the regime of a half-filled band at strong repulsion such that the energy of 
the |-|-) and |— ) states is lower than the energy of the |0) and |2) states by U/2. Under 
strong repulsion, it is reasonable to proceed from the so-called spinless superalgebra, which 
is constructed as follows. 

1. We start with the superalgebra of Hubbard operators, denoted by 5*2. From the 
eight Fermi-like Hubbard operators, we form creation and annihilation operators taking two 
orientations of the spin a into account: 

a-,. = X°+ + X~2, a+ = X+° + X2-, 

(10) 

a^ = X°--X+2, a+ = X-°-X2-. 

Their anticommutation relations are 

{aj:, 0^'}+ = a^aa' + aa'a^ = 5„a'- (11) 

These operators differ from the Hubbard operators in their (anti)commutation relations and 
therefore give a new algebra, which together with the density and spin operators describes 
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excitations in the metallic state. In [4], this algebra was denoted hy SI. It differs from the 
algebra 5*2 of Hubbard operators. 

2. We introduce spinless operators a and a'^ as 



V2 ' 



1 



V2 



a'^a. 



(12) 



The operators (a+,a, p) constitute an algebra, denoted here by SO, because we have the 
(anti)commutation relations 



{a, a+}+ = i, [p, a+] = a+, [p, a] = -a. 



(13) 



In a matrix representation of this algebra, we construct a product of the set of dynamic Bose 
and Fermi fields with the corresponding generators. Just this product enters the nonlinear 
representation of the supergroup via the exponential map of the form 



U = exp 



X* 
X* 
\ 



X 

H- 
-X* 



X 

H+ 
X* 



\ 

-X 
X 



(14) 



where and are field variables and x X* ^^e Grassmann variables. 

3. We now form a spin vector s — {s'^, s", Sz) from the operators introduced in item 1: 



: a+at. 



V2 ^ V2 
The commutation relations for the spin operators are 



[sz,s ]_ = -s~ 



(15) 



(16) 



These operators constitute a spin algebra for the spin 1/2. 

The spin operators are given by a 2x2 matrix; in the four-dimensional square matrix 
in (14), the product of the magnetic field vector times the spin vector yields a matrix of the 
form 

/O 0\ 

Hz H+ 
H- -Hz 
\0 0/ 

4. Prom the operators introduced in item 1, we also form a vector p= (p"*", p~ , Pz) of the 
density operator: 



1 



Pz = 7r(at«t ~^t^^^- 



These operators constitute an algebra with the commutation relations 

[P^, P~]- = Pz, [Pz, P^]- = P^, [pz, P~]- = -P' 



(17) 



(18) 
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In the matrix representation, this algebra yields the matrix of density fields 











\ 


























\e- 








-Ej 



A second algebra of density operators with the operator components p — (a^-aj,, a^d^ i diCil — 
) and the same commutation relations is also possible. 
We thus formed three algebras from the Fermi-like Hubbard operators: one fermionic 
spinless and two bosonic algebras. The fermionic algebra describes creation and annihilation 
of spinless electrons at the sites, the spin algebra is used for taking the interaction of the 
magnetic moments of electrons into account, and the density bosonic algebra is designed for 
describing the interaction of density fiuctuations, for example, for the analysis of density 
charge waves. 

4.1 Deformation and contraction of a Lie algebra 

Contraction of a Lie algebra is a limit operation of a Lie algebra deformation. Because 
these techniques have not been previously used in the Hubbard model, we give the necessary 
definitions below, referring to [7]-[10]. 

Let G be the Lorentz group, g be the Lie algebra of the Lorentz group, and a be the 
subalgebra of the rotation group of the three-dimensional space. Then a contraction of g 
gives the Lie algebra of the inhomogeneous rotation group. Indeed, let {cf^} be the set of 
structure constants of the Lie algebra g m. a, fixed basis ei, . . . , e„, and let A{t) be a curve in 
the group of nondegenerate linear transformations of the g group space such that ^(1) = E. 
Let ei{t) — A{t)ei and c^j{t) be structure constants of g in the basis {ei{t)}. If the structure 
constants change as t varies, then the algebra is said to be deformed. If some of the structure 
constants vanish as t tends to zero or infinity, then such a process of changing the algebra is 
called a contraction [7]-[10]. 

4.2 Local supergroup in the Hubbard model 

The Hubbard operators depend on the index parameterizing the coordinate, i.e., they are 
given by direct products of local copies of a set of generators defined at a selected site. We 
note that this set of coordinate-independent generators constitutes a global superalgebra 
(i.e., one defined on the entire space). It can be used to construct the local supergroup in 
the Hubbard model by an exponential map. The dynamics of the system is then given by the 
dynamic Fermi and Bose fields in the exponential representation. We consider the spinless 
case in what follows. 
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5 Even Hubbard operators as generators of the group 
of four-dimensional rotations 



We introduce two triples of operators: {X°^ - X^'^, X'^^, X^O) and {X++ - X", X+-, X-+). 
The commutators of these operators are expressed only in terms of operators from the 
same triple. This means that the operators divide into two groups, which commute with 
each other, i.e., define either a direct product of two SU{2) groups or the direct product 
SU (2) (8) SU (1, 1) of two subgroups. Naturally, the choice of one of these two types of groups 
must correspond to a certain choice of the fields. 

Wc note that if the commutation relations in the second subalgebra are chosen such that 
it becomes the su{l, 1) algebra, then this leads to a factor i appearing in front of some of 
the fields in the functional integral (a Wick rotation), and the group and the space hence 
become hyperbolic. 

If linear combinations (half-sum or half-difference) of the generators of these two groups 
are taken, then it can be easily verified that their commutation relations coincide with the 
commutation relations of the four-dimensional Lorcntz group [11]. Depending on whether 
hyperbolic rotations are introduced, this rotation group can be made into the group of rota- 
tions of the four-dimensional Euclidean space. We obtain the structure of the commutation 
relations for the six generators: 

[Li, Lj] = iSijkLk, [Li, Kj] = iSijkKk, [Ki,Kj] — ieij^Lk. (19) 

Here, Lj, i = 1,2, 3, specify the spatial rotation generators, Kj, j = 1, 2, 3, are the generators 
of the time axis rotations, and £yfe is the Levi-Civita tensor. Each of these groups constructed 
on the given generators is a three-parameter group. Three parameters, multiplication by the 
rotation generators in the exponential map, are responsible for the rotations, and three are 
responsible for the hyperbolic rotations (boosts). Hubbard operators are proportional to 
the following combinations: the spin operators Si = Li + Ki, and the density operators 
Pi^ Li - Ki, i = 1, 2, 3. 

We hence conclude that the Bose Hubbard operators can be divided into two subgroups, 
from which the full set of the Lorentz group generators in a spinor basis can be constructed. 

Wc note that the Lorentz group is a subgroup of the conformal group, which is the 
largest nonlinear group of the four- dimensional space-time on which the wave function of 
the system is defined. Thus extending the group of rotations of the four-dimensional space to 
the conformal and further to the superconformal group, we can construct a local (super)group 
of the Hubbard model. Because the atomic basis is a central point in practically all of the 
strongly interacting systems, we have revealed a hidden symmetry in strongly correlated 
system and indicated a way to introduce local (super)groups for the popular models. 

Another important point refers to the introduction of a matrix representation for shift 
generators. The shift generators are usually given by differential operators (derivatives). 
In the Hubbard model, just a matrix representation is necessary, different from the one 
that follows from the affinization procedure [5]. We introduce shift operators in terms of 
decreasing the dimension of the rotation group and the procedure of deforming the algebra 
of the Lorentz group generators in the four-dimensional space. This procedure is more 
complicated than affinization, but we believe it is the only one applicable to the Hubbard 
model. 
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We introduce a deformation procedure and then a contraction of the even Hubbard 
generators. We multiply the generators of hyperbohc rotations by a constant R, K'- — RKi, 
i — 1, 2, 3, such that the commutation relations become 

[KK;^]^t(^-^y.,kL,. (20) 

In what follows, similarly to [12], we regard Lj as the generators of rotations in the three- 
dimensional space and K'- as the generators of translations. It follows from (20) that via such 
a transition from the four-dimensional to the three-dimensional space, we have effected a 
transition from the group of four-dimensional rotations to the group of inhomogeneous three- 
dimensional rotations, which contain translations in a curved space. In the limit as i? — )■ oo, 
we obtain the three-dimensional Poincare group I SO (3), with the algebra of generators 

[Li, Lj] — ieijkLk, [Li, Pj] — ieijkPk, [Pi, Pj] — 0, (21) 

where Pi = Kj are generators of translations in flat space. The boost generators thus define 
translations in a curved space whose radius is defined by the expression E"^ — E^+E'^E~, and 
their identification with the generators of translations in the curved space allows introducing 
the coordinates of the base space by replacing [E"^ , E~ , E^) with {x,y,z). We see in what 
follows that this expression naturally enters the matrix elements of the supergroup and 
determines the radius of the space where magnetic dynamic fluctuations occur. 



6 The group space 

The group structure that we have revealed on the even Hubbard generators has a subgroup 
and a quotient group. Based on this, we can introduce the group of motions specifying the co- 
ordinate functions and the functions that are dynamic fields in our problem [3]. We consider 
the finite continuous group G depending on the parameters E^_. E^, E~ and Hz, , H~ . In 
what follows, we use the notation E^, E^ , E~ with the understanding that these can be re- 
placed with {x,y,z) whenever necessary. The parameters {a,h} = {E^, E'^, E~; H^, H'^, H~} 
can be considered the coordinates of a point in the (3-|-3)-dimensional group space. Each 
point of the space is set in correspondence with a transformation from the group G. The 
point corresponding to the identity transformation is called the initial point of the space. 
The initial and an arbitrary point of the space define a vector. Any infinitesimal transfor- 
mation is expressed in terms of the generators of the quotient group and the subgroup 

dG ^i{da''Xk + db"Ya), (22) 

where a'' and 6" are parameters of the group, k = 1,2,3, and a = 1,2,3. We again note 
that the quotient group here specifics the "true" coordinates {x,y,z). 
The 1-forms are expressed in terms of the group elements as 

G dG-^ = G-^ dG = i{uj'Xi + ry^), (23) 

where Xi is a generator of the quotient group [13]. 
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7 Structure equations for the group space 



Structure equations for the group space coincide in form with the Maurer-Cartan structure 
equations for a Riemannian space with zero torsion and a nonzero curvature [14]-[17]: 



(24) 



The forms cu* are components of an infinitesimal shift of the origin of a frame with respect 
to the frame at a point a. The forms cu* are a change in the components of the frame 
itself. A transformation from the group G is a rotation if it belongs to the subgroup H. The 
transformations from the subgroup H leave the origin of the group space fixed and constitute 
the so-called stationary subgroup of this space. A transformation is a shift if it is generated 
by an infinitesimal transformation cu'^Xi. 

A general transformation of the group G is represented as the product 

G^K{a)H{b), (25) 

where K{a) = e*"^"^J is a transformation belonging to the left cosct class G/H of the group 
G with respect to the subgroup H and H{b) = e*''"^" is a transformation belonging to the 
subgroup H. 

8 Cartan forms for the Bose Hubbard operators 

We consider vectors (x, y, z) of the coordinates of the S0{3) or S0(2, 1) group space, which 
in matrix (14) are expressed as components of the electric field Ez, E^, E^ , the magnetic field 
Hz{x,y,z), H~^{x,y, z), H~{x,y, z), and the Grassmann x(a;,|/,2;) and conjugate x*(a;, j/,2;) 
fermionic fields as functions of the space-time coordinates (x, y, z). We use the local spherical 
coordinates {E, 9, cp) and (if, cr, i?) to write the electric and magnetic field vectors: 

E, = Ecos9, E+ = EsmO ■e''^, = Ssin^ • e"*'^, 

H.^Hcosa, H+ ^ Hsma -e'^, ^ H siua ■ e''^ . 

Then 

E'^^El + E^E-, ^ H^ + H-^H-. (26) 

We place these fields into a supermatrix and write the exponential map in the form 



U — exp 



^Ez xi X2 E+\ 
Xl Hz H+ X3 
X*2 H- -Hz X4 



(27) 



Here and hereafter, in the expression for the exponential map, we do not write the standard 
minus sign in front of the matrix in the exponent. This sign can easily be compensated at 
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the final stage by redefining the fields. As a result, we obtain a nonlinear representation for 
the supergroup that locally performs superrotations of the atomic basis. The problem of 
calculating the matrix elements of this supergroup arises. We note that this problem is very 
complicated in the general case. In [18], we could evaluate the chiral representation, when 
matrix (27) involves only the Grassmann fields marked with an asterisk and bosonic fields. 
Below, we keep to the following strategy. 

We first evaluate the matrix consisting of only the components of the electric field E and 
the magnetic field H, whose contribution to the effective functional is zero in tlie fermionic 
fields: 



exp 





H+ 



(E, 

H 

H- 

\E- 

, ^ ^ sinhE' 
cosh E + E. 





E 



E+\ 



-EJ 



cosh H + H, 



sinhiJ 



E- 





sinh E 
E 



H 



sinh H 
H 



H 





sinh H 



E+ 



sinhE' 
E 



\ 



H 



cosh H — H. 



sinh H 
H 






cosh E — E. 



sinh E 



E J 

(28) 



The Cartan differential form is given by 






^41 





1^22 ^23 

^32 1^33 





Wi4\ 




a;44/ 



(29) 



or 



+ L02J,^~^ + <^32^ + ^^33^ 



(30) 



We consider a three-dimensional curved space with the curvature radius defined in terms of 

the invariant E. A similar argument can also be presented for the H vector, which means 
their independence of the coordinates; we therefore assume that E = const, dE = 0, and 
H = const, dH = 0. This leads to the constraints for the components of the E and H 
vectors 



dE'^ = 2E, dE, + E+dE- + E-dE+ = 0, dH^ = 0, 



(31) 



once again indicating the curved nature of both the base space and the space where the 
dynamic fields exist. 
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For convenience, we represent the Cartan 1-forms as matrix elements of two 2x2 matrices: 
we have 







\OJ32 


^^33 J 



sinh H , sinh^ H , , sinh^ H , , sinh H \ 

gi^^dH, ——H+dH- -—^H+dH, + gi^^dH+ 

sinh H , sinh H sinh H , sinh^ H 



\ j^H-dH, + g,^^dH+ -g^——dH, - -^^H-dH+J 



(32) 



where 



o-H^ = — — dx + -Tz—dy + -7:—dz 
ax ay oz 



and similarly for dH^, g2 = coshi/ — sinh H/H, and gi = coshiJ + Hz sinh H/H. The 
matrix ( ) for the electric field components is entirely equivalent to matrix (32) with 
the replacement — > E^, — > E^. In what follows, we use expression (32) to obtain the 
effective functional for the dynamic magnetic fields. 

The Cartan forms for the fermionic Hubbard generators are formed similarly to the 
bosonic case, with the only difference that odd Grassmann fields and the corresponding 
differentials are involved and the inverse matrix is calculated differently. In the final analysis, 
these Cartan 1-forms yield a functional of the fermionic dynamic fields. 



9 Effective functional 



The construction of a "symmetry tower" was proposed in [4]. It is based on the following 
observation: the variables in strongly correlated systems, being generators of some superal- 
gebras under a variation of the Hubbard repulsion and the chemical potential, constitute a 
chain of superalgebras (5*2, 5*1, 5*0) that specify excitations of those phases that emerge as 
the repulsion increases. This chain, called the "symmetry tower," is characterized by the 
number of odd generators of the superalgebras and starts with a minimal superalgebra that 
has two fermionic generators of the algebra of spinless fermions. 

In this section, we proceed from this superalgebra and construct its deformation, which is 
nonlinear in the (super) generators involved. We then evaluate the nonlinear exponential rep- 
resentation of the obtained superalgebra. Next, using formula (1) with Hamiltonian (3), we 
obtain the effective functional for the fields taking values in the chosen supergroup. We note 
that the quantum super- Yangian first occurred as an example of a nonlinear superalgebra in 
the one-dimensional Hubbard model in [19]. 

We take the expression for the supergroup yielding the supercoherent states \G) — U\0) 
for the specified regime observed in the Hubbard model under strong repulsion and half-filling 
in the form 

X X 
Hz 

H- -Hz 



U — exp 



(Ez 

r 

X* 



E+ 
X 
X 

-Ez/ 



(33) 



The generators (c+, c, c'^c, and 1 — 75/2) involved here are equal to the corresponding 
matrices at the dynamic fields in the exponent and form a closed superalgebra of a nonlinear 
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type. We give matrix expressions for the relevant generators: 



/O 0\ 

10 
10 

\o 1 1 oy 

/o o\ 
Olio 

Olio 
VO 2/ 

The (anti) commutation relations are 

{c+,c}+ = 2(l + s+ + 0, 
^+ 1 - 75 1 - 



c^c 



/O 1 1 0\ 

1 

1 

\0 0/ 

/O 0\ 



1-75 



10 

10 

yo oy 



c , s ' 



c ' c, s 



2 2 

, 1 - 75 1 + 



2 2 

^ 1 - 75 1 - 5^ _ 
2 2 

1 - 75 1 + 



_ 1 - 75 1 + ,s, 
2 2 

_ 1 - 75 1 - 
2 2 

1 - 75 1 + S;j 



2 2 

1 - 75 1 - 



2 



2 2 
[c+c,s-]_ 



-c, 



-c, 



+1 1-75 1-75 
' P I = c— — c, 



[c+,p-]=0, [c+,p,] = c^ 



c, P ] =0, [c, p ] = c 



1 + 75 1 + 75 

—2 2-' 

[c+c,p-] = -2p-, 



[c, Pz] = -c, 
[c+c, p^] = 0. 



(34) 



(35) 



We do not take the coefficients of the matrix representation of operators (15) and (17) into 

account here. 

Calculating a nonlinear representation of the above (supcr)group is cumbersome and 
quite laborious. The problem of calculating the inverse (super) matrix, which is required for 
calculating the Cartan (super)differential forms, is especially complicated. In the calcula- 
tions, we used a computer with 2 Gbytes of RAM and the Mathematica 7.0 symbohc calcula- 
tion system extended by the superEDC program [20] for working with Grassmann numbers, 
(super)matrices, and Cartan (super)differential forms. The (super)matrix elements for the 
matrix U in (33) and the Cartan 1-forms were evaluated in [21], where all the notation used 
here was also given. With the stabihzed algorithm for working with superfields that take 
values in (super)matrices, the total time for computing the effective functional given below 
was about 50 hours. 
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In calculating functional (1), we used the expansion of the kinetic energy in formula (3) 
and a representation of the product of operators at the neighboring sites in the form of a 
series in the lattice constant. In the first order, we obtain 



aUWU + b) = a{j) exp(^b^^a+{j) 



= a{x)a^{x) + a{x){b ■ V)a+(a;) + . . . . (36) 

Wc replace the index j with the coordinate x = bj for future use in deriving the repre- 
sentation of the effective functional. We note that the operators a and a"*" are expressed in 
terms of the Hubbard generators. 

Expanding the exponential in a series using formula (36) and the expression for the 
Cartan (super) differentials in the first order in the lattice constant, we obtain the effective 
fermionic functional in the form of the sum 

Li = ih + h)x{r){b ■ V)x*(r) + hx{r)x*{r). (37) 

This functional is quadratic in the spinless fermionic field and contains the dynamic fields 
Hi and a coordinate dependence that gives the metric of the curved space. The coefficients 
ki in (37) are given by 

ko^EH[f^ + f2hs + f2{x + z)], 

ki^f2 + fk + V'lhhihi + x + y + z) + fUHzh2 + zh2 + {x + y){H, + z)], 

(38) 

k2 = ko[E^{h - f2h,) + (y + z){zh - E^f2) + {H, - z){U - fsh^)], 

k3 = ko{2{h - h,f) + {x-y- h^){h - Kf2) - {v + z)[2f, + {x - y)f2]}. 

These expression include contributions from the derivatives of the magnetic fields; the exact 
formula for them will be given elsewhere. The other elements of the condensed notation are 
presented in the appendix. 

Matrix (32) gives differentials of the dynamic fields. We can see from (1) which derivatives 
of these fields appear. It follows from (1) that the effective functional involves only the 
expression for 0022- With all the foregoing, we obtain the expression for the effective functional 
of bosonic fields: 

Lo^-[H-{b-VH+) + H+{b-VH-) + H,{b-VH,)]^^^^^-^ + H,{b-VH,f-^^^c^^^ (39) 

If the second-order terms in the lattice constant are kept in the kinetic energy expansion 
in (36), then an effective functional can be obtained that contains quadratic contributions 
from derivatives of the magnetic fields. This then gives a nonlinear sigma-model, which is 
an analogue of the Heisenberg model in the functional integral for the Hubbard model. A 
functional linear in the derivatives is similar to the functional in the Chern-Simons model, 
which is natural for our problem. 
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Finally, the full effective functional for the spinless version of the Hubbard model has the 
form 

Leff = 4 + ^l. (40) 

We believe that the regime we have chosen in the Hubbard model is related to the spin- 
liquid phase revealed in the organic metal (BEDT-TTF)2Cu2(CN)3. This compound requires 
taking a two-dimensional lattice with a hexagonal symmetry. 

10 Conclusions 

We have offered a further development of the approach formulated in [3], [4]. We showed 
that the supercoherent state of form (33) defines a nonlinear representation of the minimal 
(super) extension of the Lorentz group. Continuing this procedure of (super)extension leads 
to a representation of the superconformal group in the superspinor basis, nonlinear in the 
fermionic fields. We constructed a new nonlinear matrix representation for the translation 
generators in the (2+l)-dimensional space-time. It defines a nonlinear super algebra that 
enlarges the set of (super) algebras that have been proposed in the "symmetry tower." We 
note that after [19], the search for a quantum algebra in the Hubbard model was conducted 
very vigorously but unfortunately unsuccessfully. We have indicated a way to introduce 
similar structures in the Hubbard model. Superalgebra (35) contains pairs of odd generators 
that anticommute with the shift operator (for a spherical space, a shift is equivalent to a 
rotation), as is appropriate for the standard supersymmetry algebra. 

The effective functional of form (40) explicitly contains the coordinates of the S0{2, 1) 
or 5*0(3) group spaces and describes the dynamics of Bose and Fermi fields in a three- 
dimensional curved space-time. A functional of this type, as far as we know, was first 
obtained in the Hubbard model. 

In relation to this work, we note an entirely new problem for the Hubbard model: cal- 
culating the cohomology groups, the second and the third in particular, with the purpose 
of studying the superalgebras in the "symmetry tower." This approach opens the possibil- 
ity of studying the mechanism of the Mott-Hubbard-type metal-insulator transition as a 
spontaneous violation of the (super) conformal (super)group. 
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Appendix 

The additional notation in formulas (38) is 

{sinh{aE))/E-{sinh{aH))/H _ d , 

J[a) E'^-H^ ' ~ d(lP~y /(«)J«=i> 

d d 

flhh = -Q^{fhh)a=l, f2h = 'gff^^f'^^^ 

(sinh E)/E — (sinh if) /H cosh E — cosh H E sinh E — H sinh H 

^2 — _ JJ2 ' ~ J^2 _ JJ2 ' •^^ ~ J^2 _ JJ2 ' 

, _ sinh H fi — ^^"^^ ^ f f — ^^^^ ^ f f — ^^^^ ^ f 

Jhh ~ 'J^ fhh, flhh ~ flhh, 72 — 'J^ /2; /2/1 — /2h; 



ho = H+-H-, hi = H+ + H-, h2 = H~ + H, 
hs = H — Hz, /i4 = + Hz, = H^ — H. 



Z1 
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